Propagating quantum coherence for a biological advantage 
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Photosynthetic antennas need to efficiently transfer absorbed energy through chains of pigment molecules 
to reaction centers, where excitations induce charge separation. Experimental evidence shows that a variety 
of photosynthetic systems can preserve quantum beats in this energy transfer process, even at room tempera- 
ture. However, whether this quantum coherence arises in vivo and whether it has any biological function have 
remained unclear. Here we present a theoretical model which suggests that the creation and recreation of co- 
herence under natural conditions is ubiquitous. Our model allows us to demonstrate a quantum coherent ratchet 
effect, in a design based on an energy transfer pathway in the Fenna-Matthews-Olson complex of the green 
sulfur bacteria. These results suggest a biological role for coherent oscillations in spatially directing energy 
transfer. 



I. INTRODUCTION 

Mounting experimental evidence for electronic quan- 
tum coherence in photosynthetic energy transfer |[]}[7) has 
spawned much debate about both the detailed nature and the 
biological role of such quantum dynamical features. Quan- 
tum coherence is usually encountered in the first, light har- 
vesting stage of photosynthesis. It includes two distinct phe- 
nomena that can be differentiated by the choice of basis used 
to describe the electronic excitations. In the site basis, cor- 
responding to the excitation of individual pigment molecules, 
coherence emerges in molecular aggregates even in thermal 
equilibrium, since eigenstates are delocalized over multiple 
chromophores. There is strong evidence that this coherence, 
which under a single excitation is equivalent to quantum en- 
tanglement [8,9], can enhance biological energy transfer ifTQl — 
H3l . In contrast, it is coherence in the exciton basis, that is, su- 
perpositions of the energy eigenstates, which drives quantum 
beating via the Schrodinger equation. Photosynthetic systems 
at ambient temperatures have been shown to exhibit this kind 
of quantum beating when artificially excited fl6j 0, but the 
significance of these discoveries for biological energy transfer 
remains unclear. 

The Fenna-Matthews-Olson (FMO) complex of green sul- 
fur bacteria is an extensively studied protein-pigment com- 
plex lfl4l -tT6 l that exhibits such quantum beating at room tem- 
perature [6|. Biologically, FMO acts as an energy transmit- 
ting wire, delivering an electronic excitation created by pho- 
ton absorption in the chlorosomone antenna to a reaction cen- 
ter where it induces charge separation. These energy transfer 
pathways are illustrated in Figure[T] The path section from site 
1 to 2 is particularly unusual, since it is energetically uphill 
while these sites also have the strongest electronic coupling 
of any pair of sites in the complex. It has been speculated that 
these factors may indicate a role for quantum coherence in 
contributing to unidirectional energy flow through this system 
by avoiding trapping in local minima of the energy landscape 
ifTTl . All other steps of the FMO electronic energy transfer 
pathways in the direction towards the reaction center are en- 
ergetically downhill, consistent with a pathway that would be 




FIG. 1: Energy transfer pathways in a monomer of the FMO com- 
plex of Chlorobaculum tepidum (C. tepidum). (a) Side view of a 
monomer of the FMO complex [16|, showing the primary energy 
transfer pathways 1 1 8 1 and the inter-complex coupling (ICC) basis 
states that couple site 8 to the remainder of the complex. Site occu- 
pation probabilities for the ICC basis states are proportional to the 
area of the colored circles, (b) Site energies along the upper en- 
ergy transfer pathway depicted in panel (a), with the energy of site 8 
approximated by the antenna baseplate energy 1171 . Lines between 
sites indicate weak (dashed, 10 cm -1 < J < 40 cm -1 ) and strong 
(solid, J > 40 cm -1 ) electronic couplings. 



optimal for classical energy transport. 

This paper addresses two open questions concerning the 
relevance of excitonic coherence to photosynthetic energy 
transfer. First, we present a theoretical framework which sug- 
gests how coherence should arise and recur in the process of 
energy transfer. This is the 'propagation' of quantum coher- 
ence, whereby a process of continual renewal following inco- 
herent quantum jumps may allow non-zero coherence to last 
indefinitely, despite rapid decay after each jump. Second, we 
present evidence that such coherence is relevant to the func- 
tion of photosynthetic systems. Inspired by the uphill ener- 
getic step in FMO, we construct an explicitly solvable model 
to show that the non-equilibrium nature of even limited quan- 
tum beating may allow for qualitatively new types of dynam- 
ics. Under biologically plausible conditions, these dynami- 
cal features allow for the operation of a coherently enabled 
"quantum ratchet" to enhance directed energy flow through 
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the FMO complex. 



H. PROPAGATING QUANTUM COHERENCE 

In photosynthetic energy transfer, excitons typically need 
to travel through a series of protein-pigment complexes be- 
fore reaching reaction centers [15|. To accurately understand 
the role of coherent dynamics in individual complexes such as 
a monomer of FMO, it is first necessary to understand which 
particular states donate or accept excitations for inter-complex 
transfer. Whether or not coherence arises under natural condi- 
tions |fT9ll2"0l is informed by the nature of these states. Since 
inter-complex couplings are relatively weak, in our analysis 
we treat them perturbatively, as in multichromophoric gener- 
alizations of Forster theory used to calculate overall transi- 
tion rates between donor and acceptor complexes consisting 
of multiple chromophores lflOl4T2ll . 

Our starting point is the equation of motion for the reduced 
density matrix, which is derived with the following adapta- 
tion of the multichromophoric energy transfer rate model lfl2l . 
The zeroth order Hamiltonian is Hq = Hd + Ha where 
H D = H e D + £ tf B Di} |A)(£>il + H g D , with H% the elec- 
tronic Hamiltonian of the donor complex, and correspond- 
ing definitions for the acceptor complex A. States \Dj) and 
\Ak) for j — 1, . . . , n and k = 1, . . . , m form an arbitrary 
orthonormal basis for donor and acceptor single-excitation 
electronic states and Bp... are bath operators that couple the 
electronic chromophore states to environmental states of the 
pigment-chromophore system. The ground state donor (ac- 
ceptor) bath Hamiltonian H 9 D can be taken without loss of 
generality to be a set of independent harmonic oscillators. We 
assume that no bath modes are coupled to both the donor and 
acceptor so that [He,, Ha] = |12|. The donor and acceptor 
complexes are coupled by a dipolar interaction H c = J + 
with J = J2jk Jjk\Dj)(Ak\- Calculating the evolution of the 
reduced system density matrix a = Tib P to second order in 
H c (" iPTZl and Appendix [A| yields 
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for the reduced acceptor and donor density matrix elements, 
respectively, where E£>(£, oS) and I A (w) denote matrices of 



donor and acceptor lineshape functions (see Eqs. (A4 - A5| ). 

Instead of focusing on the multichromophoric energy trans- 
fer rate between complexes that results from summing Eq. ([T} 
or (|2) over all diagonal terms [ 12 1, we focus here on impor- 
tant features relevant to quantum coherence apparent from 
Eqs. ([TJj2]l directly. These equations indicate that accep- 
tor populations |y4fc)(y4fc| will grow and donor populations 
\Dj)(Dj | will decay only if there is at least one non-zero cou- 
pling term Jj k — (Dj\ J\A k ) to those states. Accordingly, we 



argue that the transfer of electronic states is most sensibly de- 
scribed by the "inter-complex coupling basis" basis in which 
J is diagonal, rather than by either the site or exciton (energy) 
basis. This inter-complex coupling (ICC) basis is given by the 
singular value decomposition J = Up 3U\, where J is a rect- 
angular diagonal matrix and Ud and Ua are unitary transfor- 
mations of donor and acceptor electronic states. We can thus 
write the inter-complex coupling as H c = J^i Ji(\Di) (Ai \ + 
\Ai)(Di\) in terms of the ICC states \D t ) = U D \D t ) and 
\A~i) = Ua\Ai) fox I e {1, . . . , min(n, to)}. In the ICC basis, 
the full electronic Hamiltonian in block-matrix form is 
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Since in general the transformation that diagonalizes J will 
not coincide with the (exciton) eigenbases of and H A , 
population growth of an acceptor ICC state \A{) (Ai | thus cor- 
responds to growth of excitonic coherences. 

Although in principle Eqs. ([T]-[2} specify all dynamics rele- 
vant to inter-complex transfer, the time-dependent donor line- 
shape Ed(4,w) obscures the specific nature of the depen- 
dence on the donor. Accordingly, we also derive a time- 
convolutionless quantum master equation (Appendix [Bj un- 
der the usual assumption of weak coupling to the bath |21|. 
Under this approximation, we see that growth of an acceptor 
population \A{)(Ai \ is proportional to populations only of the 
coupled donor \D{)(Di\. Likewise, decay of a donor popula- 
tion \D[)(Di | is proportional to populations only of that donor 
itself (see Eqs. (B7 B8 i in Appendix IE). Accordingly, inter- 



complex transfer rates may show oscillations reflecting donor 
quantum beats, since the ICC states on the donor which trans- 
mit excitations do not necessarily correspond to energy eigen- 
states. While this argument is only rigorous in the case of 
weak coupling to the bath, which is not necessarily the case 
for FMO and other light harvesting systems jTT), our simu- 
lations find excellent agreement even for moderate strength 
environmental coupling, as we show below. 

To test this analysis of inter-complex energy transfer, we 
first consider its predictions for the special case in which there 
is only one non-zero inter-complex coupling in the ICC basis, 
H c = JJD*)(A*\ + h.c. If the acceptor is always initial- 
ized in the state p* A = \A*)(A*\, then when back-transfer to 
the donor is neglected as is valid in the perturbative limit, the 
acceptor density matrix should be well described by 
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where Q(t) is the Greens function denoting evolution of the 
acceptor-bath system for time t, with the bath initialized at 
equilibrium. If the predicted donor state is correct, then ne- 
glecting temporary bath reorganization effects, the rate of 
inter-complex transfer should then be proportional to the pop- 
ulation of the predicted donor state, for a predicted inter- 
complex transfer rate 
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FIG. 2: Testing the theory of propagation of coherence via the inter- 
complex coupling (ICC) basis, (a) Simulated (solid black, 2CTNL) 
and predicted (dashed blue, Eq. {5}) inter-complex transfer rate for 
the coupled dimer model described in the text, with ICC donor 
D*} = |2) and the initial condition \^o) = (b-d) Population of 
site 3 (b), real (c) and imaginary (d) parts of the 3-4 site coherence 
for the simulated (solid black) and predicted (dashed blue, Eq. (4}) 
acceptor state normalized to unity, with ICC acceptor p* A = j3)(3j. 



where p^* (t) denotes the probability of the donor being in the 
state \D*) and pa the total probability of the excitation being 
on the acceptor. For a model system, these predictions show 
remarkable agreement with results derived from an indepen- 
dent simulation based on a 2nd-order cumulant time-nonlocal 
(2CTNL) quantum master equation [22 1 . We consider transfer 
between two dimer complexes (labelled sites 1, 2 and sites 3, 
4), with intra-dimer Hamiltonian parameters matching those 
of the 1-2 dimer of FMO (Appendix |EJ Eq. (JETJ) and inter- 
complex coupling J = J*|2)(3|. We perform calculations 
in the limit J* — > (see Appendix [P|, to ensure accuracy 
of the perturbative description and eliminate back-transfer ef- 
fects. The 2CTNL calculations are carried out at 300 K with a 
bath modeled by a Debye spectral density with reorganization 
energy 35 cm -1 and correlation time 50 fs ifTTl . with the ini- 
tial condition on site 1. Figure [2] compares simulated 2CTNL 
results with the predicted time-dependent inter-complex trans- 
fer rate, Eq. |5]), and acceptor density matrix, Eq. Q (normal- 
ized to unity for greater clarity), calculated from the 2CTNL 
results. The results show that estimates based on the dominant 
elements of the ICC basis provide an accurate representation 
of both the energy transfer rate (panel a) and acceptor density 
matrix (panels b-d). Transfer of excitation in the ICC basis 
from \D*) to \ A*) produces a superposition of acceptor eigen- 
states (of Ha) that gives rise to excitonic beating and hence 
to oscillatory behavior of the site populations (panel b) and 
coherences (panels c-d). Two features are of particular signif- 
icance for inter-complex transport. First, that the inter-dimer 
transfer rate clearly tracks coherent oscillations of the donor 



population | D* ) (D* | (panel a) and second, that the acceptor is 
initialized in a state with non-zero excitonic coherence (panels 
b-d). 

A simple example of the usefulness of the ICC basis is to 
determine the initial conditions for electronic excitation trans- 
fer through the FMO complex. The recently discovered 8th 
chromophore ||23l l24l provides a plausible donor to the re- 
mainder of the complex since it sits on the side nearest the 
chlorosome antenna complex [25]. Since structural informa- 
tion concerning the location of the FMO complex is limited, 
standard practice to date has been to choose initial and final 
conditions for simulation of energy transfer in FMO based on 
approximate orientation and proximity of chromophores. The 
choice of such initial conditions has varied [14, 25, 26], partic- 
ularly with regards to whether or not the initial quantum states 
include any excitonic coherence. Evaluation of the ICC basis 
between a donor complex consisting solely of site 8 and an 
acceptor complex consisting of the remainder of the complex 
(sites 1-7) implies that the acceptor state is mostly localized on 
site 1 (see Appendix IB. We illustrate this in Figure [Tf a). This 
initial condition is not an energy eigenstate, so the resulting 
in vivo dynamics would necessarily give rise to the quantum 
beating seen in the laboratory experiments 0. 



III. COHERENT VERSUS THERMAL TRANSPORT IN A 
MODEL DIMER WITH AN ENERGY GRADIENT 

The way in which this coherence regenerating transfer 
mechanism can yield a non-trivial biological role for coher- 
ence can be already illustrated with a simple model dimer 
complex connected to other complexes. Choosing inter- 
complex couplings to and from the dimer to be at individual 
sites as above implies that initialization and transfer in the ICC 
basis will be at these sites. Consider preferred 'forward' exci- 
tation transfer to be that from the donor at site 2 onward to the 
next complex, so that the asymptotic probability of success- 
ful transfer through the complex will be proportional to that 
population. For a dimer, the electronic Hamiltonian is given 
by 
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where 9 is the mixing angle, which measures the intra-dimer 
derealization, and AE the exciton energy difference. A 
non-zero mixing angle corresponds to non-zero exciton de- 
localization, as indicated by the inverse participation ratio 
N = l/(sin 4 9 + cos 4 9). 

The dimer admits two extreme models of energy transfer: 
quantum beating due to coherent evolution and instantaneous 
relaxation to thermal equilibrium between excited electronic 
states. Instantaneous relaxation provides an upper bound on 
the speed of excitation transfer governed by a classical mas- 
ter equation, since the dynamics governed by such equations 
are driven toward thermal equilibrium. This is imposed by 
the requirement of detailed balance which governs classical 
dynamics, whether between sites as in Forster theory ll27l or 
between exciton populations as in variants of Redfield theory 
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FIG. 3: Difference between thermal and coherent population on site 
2, P2 ~ as a function of dimer Hamiltonian parameters for (a) 

initial site i = 1 and (b) initial site i — 2. The empty symbols □ 
and o indicate location of parameters for the 1-2 dimer in the FMO 
complex of C. tepidum at 300 K as determined in Refs. H4landl3T1 
respectively. Filled symbols indicate the corresponding parameters 
at77K. 



For instantaneous relaxation to the thermal distribution, 
the probability that site 2 is occupied is independent of initial 
condition: 



1 



In contrast, boosts in population due to quantum beating are 
not restricted by such classical limits 1221 |29l |30l . For co- 
herent motion with initialization at site 1, the time-averaged 
probability of an excitation at site 2 is 



Pi -^2 = (|{l|e-^ t |2>r) f = 2cos 2 



' sin 



(7) 



while for initialization at site 2, we have — 1 — Pi° t >2- 

Figure [3] plots the difference between coherent and thermal 
population on site 2, as a function of both the intra-dimer de- 
localization measure 8 and energy difference AE, for both 
possible initial conditions. It is evident that regardless of ini- 
tial conditions, for a sufficiently uphill energetic arrangement 
(AE > 0) intra-dimer quantum beating will be asymptot- 
ically more effective than intra-dimer thermalization in en- 
abling transfer onward from the complex via site 2. The loca- 
tion of the FMO parameters in Fig.[3]shows that the 1-2 dimer 
of FMO satisfies such an arrangement at 77 K and is on the 
borderline for coherence enhanced transfer at room tempera- 
ture. 



IV. DESIGN FOR A QUANTUM COHERENT 
BIOLOGICAL RATCHET 



The asymmetry between incoherent and coherent popula- 
tion transfer seen above for a simple model dimer suggests a 
design principle that could be exploited for enhanced unidi- 
rectional transfer [ 17 1 and, more generally, as a novel type of 
quantum ratchet 1 32 1 . Ratchets and Brownian motors 11331 uti- 
lize a combination of thermal and unbiased non-equilibrium 
motion to drive directed transport in the presence of broken 
symmetry. To take advantage of such a ratchet effect, strongly 
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FIG. 4: Biased energy transport in an excitonic wire due to propagat- 
ing quantum coherence, (a) A weakly linked chain of heterodimers is 
arranged such that the higher energy state is always to the right, with 
a typical inter-dimer distance of 3 nm. The arrow indicates the di- 
rection of biased transport, (b) Relative asymmetry between left and 
right inter-dimer transfer rates (Eq. {D19]l of Appendix[D| as a func- 
tion of the time before transferring for a dimer excitation initialized 
in the asymptotic distribution of site populations, (c) Drift velocity 
vs coherence time as modified by bath correlation time (squares) and 
cross correlation coefficient between dimer sites (diamonds). The 
dashed line is a linear fit to guide the eye. Full details of the simula- 
tions in panels b and c are in Appendix [D] 



linked chromophores with coherent transfer not limited by de- 
tailed balance should have an uphill energy step relative to the 
desired direction of transport, whereas weakly linked chro- 
mophores with incoherent transfer steps should be arranged 
downhill. 

As a proof of principle, we present a dramatic example 
in which this coherent ratcheting effect results in asymptotic 
spatial bias of transport. Consider a weakly linked chain of 
heterodimers breaking spatial inversion symmetry, as illus- 
trated in Figure |4^ a )- I n anv classical random walk, transi- 
tion rates must satisfy detailed balance to assure thermal equi- 
librium, which guarantees that a classical walk along such a 
chain is unbiased (Appendix [CJ. However, we have carried 
out 2CTNL quantum simulations on small chains of dimers 
that suggest including the effects of coherence in each dimer 
breaks the symmetry of detailed balance to yield a non-zero 
drift velocity. Since simulations with the 2CTNL approach 
would be computationally prohibitive for large numbers of 
dimers, these simulations were carried out on a chain of three 
weakly linked dimers with parameters for each dimer match- 
ing those of the 1-2 dimer of FMO used earlier and an inter- 
dimer coupling of 15 cm -1 . We note that this inter-complex 
coupling strength is well below the cut-off below which en- 
ergy transport in light harvesting complexes is usually de- 
scribed by completely incoherent hopping within Forster the- 
ory, though without the possibility of coherence regeneration 
lfT8l [3D . The results of these simulations are used to de- 
fine left and right inter-dimer transition rates for the central 
dimer. These are then used together with ICC initial condi- 
tions from our analysis of the inter-complex coupling as input 
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into a generalized classical random walk describing energy 
transfer along a chain of dimers. The physical model corre- 
sponds to the chain shown in panel a of Figure|4] with red sites 
(dimer internal site 1 ) at energy zero, while blue sites (dimer 
internal site 2) are at energy 120 cm -1 . Formally, this theo- 
retical approach corresponds to using the microscopic quan- 
tum dynamics within and between complexes to define state- 
specific rates between complexes that generate 'quantum state 
controlled' incoherent energy transfer dynamics over long dis- 
tances. Full details of the construction of transition rates and 
of the set up and solution of the generalized random walk are 
described in Appendix [D] 

These quantum state controlled incoherent dynamics can 
generate a significant bias in the spatial distribution of excita- 
tion transfer when the intra-dimer dynamics display long last- 
ing quantum coherence. We analyzed the random walk with 
both Monte Carlo simulations on long finite chains and an an- 
alytic solution 1341 of the asymptotic mean and variance of the 
distribution, as detailed in Appendix|D] Figure^b) shows that 
the underlying asymmetry of transfer rates and violation of de- 
tailed balance dynamics is due to non-equilibrium state of the 
donor. The bias is in the forward direction, corresponding to 
the uphill step within dimers. Figure |4|c) plots the asymptotic 
drift velocity of the random walk against the timescale of ex- 
citonic coherence. We determine this coherence time from the 
timescale of exponential decay of intra-dimer excitonic coher- 
ence. The timescale of coherent oscillations was tuned in two 
ways, (i) by changing the bath correlation time and (ii) by 
increasing the spatial correlations between the chromophore- 
bath couplings |9|. We see a strikingly close correlation be- 
tween the timescale for quantum beating and the magnitude 
of the bias, regardless of the underlying physical mechanism 
used to tune the coherence time. 

These results for an infinite chain of heterodimers dramati- 
cally confirm the effectiveness of quantum coherent ratcheting 
of energy transfer from the combination of intra-complex ex- 
citonic coherence within dimers and an uphill intra-complex 
energy gradient. The non-zero drift velocity means that over 
long distances this ratchet offers a quadratic improvement in 
transfer times over any corresponding classical walk, which is 
unbiased (Appendix |C]>. However, in contrast to the speedup 
offered by quantum walks [ 35 1, this quantum coherent ratchet 
requires only short ranged and short lived coherences that will 
be resilient to the static and dynamic disorder of biological 
environments. Most importantly, this spatial bias constitutes a 
preferential direction for the energy flow across multiple light- 
harvesting complexes. It is thus of direct biological relevance 
for FMO, which serves as a quantum wire connecting the an- 
tenna complex to the reaction center. We therefore now con- 
sider the implications for FMO in more detail. 



V. RATCHETED ENERGY TRANSPORT IN THE 
FENNA-MATTHEWS-OLSON COMPLEX 

We now specifically consider the role of the coherent dy- 
namics in the uphill step energy of the FMO complex, which 
corresponds to the 1-2 dimer in the usual notation (see Fig . [TJ . 
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FIG. 5: Simulations of FMO dynamics at room temperature, (a) 
Population of dominant ICC donor state at site 2 in the 1-2 dimer 
(|Z?2}(£ ) 2|, solid line) compared with the rate of population change 
of its ICC acceptor state in the 3-7 complex (IA2XA2I, dashed line) 
and the rate of population change of the other ICC acceptor state not 
coupled to this donor state (| Ai) {^4i | , dash-dotted line), for the ini- 
tial condition is \tpo) = !!)■• The time derivatives have been scaled 
to aid comparison of correlations, (b) Population of site 2 relative 
to the total 1-2 dimer population, p2/(pi + P2), for both choices of 
initial conditions. 



Since the 1-2 dimer is relatively weakly coupled to the other 
chromophores in the complex, we may consider transfer to 
and from this dimer on the basis of our perturbation analysis 
using ICC states. By performing a singular value decompo- 
sition of the appropriate coupling matrices (see Appendix [E}, 
we find that the dominant couplings to and from this dimer are 
from site 8 to site 1, and from site 2 to site 3. This suggests 
the relevance of our dimer model from Section [IlIJ where site 
8 acts as a donor to the 1-2 dimer, and site 2 in turn acts as 
a donor to the 3-7 complex. Figure |5ja) presents results of 
a 2CTNL simulation on sites 1-7 of FMO, partitioning FMO 
between donor state on the 1-2 dimer and acceptor states on 
the remaining sites 3-7. (i.e., neglecting the prior donation 
from the 8th site to the 1-2 dimer). The corresponding ICC 
donor/acceptor states are given in Appendix|E](Eqs. ( |E5f|E6[ )). 
We see that the ICC donor population \D 2 )(D2\ ~ |2)(2| is 
positively correlated with the rate of growth of its coupled ICC 
acceptor \A 2 ) (A 2 \, but negatively correlated with the growth 
of the other ICC acceptor state | J 4i)(^4i|, to which it is not 
coupled. This is in agreement with the predictions of our the- 
ory from Section|Il] (The small deviations arise because FMO 
is not quite in the regime of validity for the perturbation theory 
and because Eq. |5]) is not strictly valid for the situation with 
two acceptor states \A).) The simulation is performed for a 
bath correlation time of 50 fs as described previously lfl7l . 

As indicated by the location of the FMO 1-2 dimer Hamil- 
tonian parameters in Figure [3] this particular chromophore 
dimer appears be optimized to have an uphill energy gradient 
just large enough so that excitonic coherence enhances trans- 
fer if initialized at site 2 (panel a) without also suppressing 
transfer initialized at site 1 (panel b). In Figure |5|b) we com- 
pare the portion of dimer population on site 2 from 2CTNL 
calculations with the classical upper bound of the thermal av- 
erage, for initial conditions in both of the ICC states |1) and 
1 2). The populations show quantum beating deriving from 
partially coherent motion. The populations also agree with 
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the predications of our simple dimer model (Section Pit: the 
population at site 2 averaged over quantum beats (Eq. (|7]i) is 
nearly equal to the thermal average (Eq. |6])) for initialization 
at site 1 and greater than the thermal average for initialization 
at site 2. 

As evident from Fig. [5] these temporary boosts in popula- 
tion at site 2 due to coherence should correspond to enhanced 
biological function, since they drive excitations preferentially 
toward site 3 instead of backwards toward the antenna com- 
plex. The 1-2 dimer in the FMO complex thus appears to act 
as one link of our proposed quantum coherent ratchet for en- 
ergy transport, thereby enhancing unidirectional energy flow 
through this pathway of the FMO complex. 



VI. CONCLUSIONS 

We have proposed a microscopic mechanism for the propa- 
gation of excitonic coherence in energy transfer between pho- 
tosynthetic complexes. Our analysis shows that the key role in 
inter-complex transfer is played by the inter-complex coupling 
(ICC) basis, rather than energy or site bases employed by all 
prior analyses. Based on ICC donor and acceptor states, we 
showed that coherence can contribute to enhanced unidirec- 
tional energy flow through an energetic pathway in the FMO 
complex. Since one-way transmission of electronic energy 
from the antenna complex to the reaction center constitutes 
the main function of the FMO protein in the light harvest- 
ing apparatus of green sulfur bacteria, this provides evidence 
for a biological role of coherence in this particular light har- 
vesting system. Direct calculation of ICC states should also 
help us evaluate the significance of long lasting coherences in 
bacterial reaction centers [3], since these systems are usually 
unlikely to absorb light directly [ 15 1 and thus might be ben- 
efitting from recurrence of coherence propagated from light 
harvesting complexes. 

More broadly, the power of propagated coherence is evi- 
dent from our example of a coherently enabled ratchet effect 
along a chain of heterodimers. This proof of principle model 
shows that even short-lived excitonic coherence can lead to 
large scale dynamics that are incompatible with any com- 
pletely classical description. Propagating coherence allows 
for coherent motion to play a role in transfer rates and pho- 
tosynthetic efficiency in light harvesting systems of arbitrary 
size. It also provides a scalable method for multi-scale mod- 
eling of large systems, since accurate fully quantum models 
need only to be used for tightly coupled sub-complexes, while 
transfer between sub-complexes may take the form of clas- 
sical hops with connections between states of different sub- 
complexes constrained by the ICC basis. We showed how 
these techniques may be applied with the demonstration of 
the quantum coherent ratcheting of energy transport along a 
heterodimer chain. Finally, these dynamics constitute a novel 
type of quantum ratchet that utilizes propagated quantum co- 
herence in place of a driving force and as such are also of 
more general interest. Thus, in addition to excitonic systems, 
we expect that a similar ratchet effect could be demonstrated 
in other experimental systems described by spin-boson Hamil- 



tonians, such as cold atoms in optical lattices ll36l[37ll . 
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Appendix A: Extending multichromphoric Forster theory 

Consider a system under zeroth order Hamiltonian Hq with 
perturbation Hamiltonian V. In the interaction picture p\{t) = 
e iH a t/h p(£j e -iH t/h ^ yon N eumann equation is 

f = -{m)Mt% (ad 

which can be formally integrated to yield, 

pi(t) = pi(0) - % -J dt'imt'),^')]. (A2) 

Inserting Eq. ( |A2| i into Eq. ( |A1[ ), keeping terms up to second 
order in V and transforming back to the Schroedinger picture 
yields a time-convolutionless equation of motion 

^ = - p J dr [V, [ e -«*>r/» 7e iffor/» p{t)]] _ (A3) 

where we approximated e~ lH ° T / h p(t — T ^ e lE -a r l h ~ p(j^ 
(valid to lowest order in V). Equations ([I||2| of Section|II]fol- 
low from direct evaluation of Eq. \A3\ with Ho and V = H c 
and with appropriate insertion of acceptor and donor lineshape 
functions, 

/oo 
dt e lut 
OO 

x ^[e^^iA^e-^/^Ak)^} (A4) 

E 3 D 3 (t,uj) =2 f dt'e- lult ' 
Jo 

x Tr B {e- lH ° t / h (D J ,\e lHDt/h p c D \D :i )} , (A5) 

where the initial condition p = p e D p g A is a product of an ex- 
cited donor state with acceptor in the ground state at thermal 
equilibrium fl2l . Note that our perturbation parameter is the 
inter-complex coupling H c rather than the usual coupling to 
the bath, so our results apply to any bath coupling strength. 
Explicit dependence upon t in the donor lineshape is removed 
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by the Markovian approximation, allowing the upper limit of 
this integral to be extended to infinity and assuming that the 
donor p e D is stationary. This would give rate expressions ex- 
tending equilibrium multichormophoric Forster theory [ 12 1. 

We note that the result in Eqs. ([TJ|2]l shows that these equa- 
tions do not necessarily conserve positivity, a feature hidden 
by the sum over states to determine an overall transfer rate 
Il2l . This is an intrinsic limitation of the perturbative ap- 
proach to inter-complex transfer. In particular, these expres- 
sions may predict the creation of non-physical acceptor co- 
herences of the form even without necessarily in- 
creasing both of the corresponding population terms \A k ) (A k \ 
and lAfc/) (Afc< |. For this reason, in determining ICC states, 
we explicitly only consider those states which will experience 
population growth or decay. 



Appendix B: Weak system-bath coupling 

Under an approximation of weak coupling between the sys- 
tem and bath degrees of freedom, the full density matrix can 
be factorized in the form p(t) = p^ q (j(t) between the equi- 
librium state of the bath = p 9 D p\ and the electronic state 



of the system a(t). Accordingly, Eq. \A3\ becomes a time- 
convolutionless master equation, 



da _ 1 

~dt ~ 



V,[K(t),a]\, 



(Bl) 



K(t)= [ dTTr B [ e - lHoT/h Ve lHoT/h p%}. (B2) 
Jo 

The Markov approximation is given by t —> 00, in which case 
we write K = lim^oo K (t). Since V and K are not in gen- 
eral equal, the Markovian expression is not in Lindblad form 
and thus does not necessarily conserve positivity, a standard 
limitation of perturbative derivations of quantum master equa- 
tions ET1 . 

For convenience, from now on we the Markov approxima- 
tion. Similar conclusions hold in the more general case. We 
then can write Eq. ( |B1| | in terms of the evolution of each den- 
sity matrix element as 



d(Ja 

dt 



ab _ 1 n 



bedded 



(B3) 



cd 



by defining Redfield-like tensor elements 

Rabcd = — ^ [SdbV ae K ec + 5 ac KdeV e b\ 

+ K ac V db + V ac K db . (B4) 

To evaluate our particular model of inter-complex transfer, it 
is useful to define an acceptor lineshape that only depends 
upon the bath state in the same form as the donor lineshape 
(Eq. dA5])), 



Tr B {e lHa o t / h (D r \e- %HDt/h \D ] )p 9 D } . (B5) 



Since K is Hermitian, to evaluate the model of Section |TT1 it 
suffices to calculate Kj k — (Dj\K\A k ). Using the cyclic 
property of the trace and inserting the donor and acceptor line- 
shapes, we find 



1 x 

K i k = J,. <L J i' k ' 
j'k> 



du£ 3 D 3 {u)I k A k {Lu) 



(B6) 



We can now evaluate the tensor elements in Eq. ( |B4| i that spec- 
ify the influence of donor density matrix elements on inter- 
complex transfer. The relevant tensor elements for the change 
of the acceptor elements due to the donor are given by, 



7? 



J jk „ 

kk'jofo ~ 4^ 



duj 



£%\u)lf"{u) 



+ J jok I k A k X") 



(B7) 



and for the change of the donor itself due to donating an exci- 
tation, 



£ 

kk'j" 



J 



j"k' 
4.7T 



dco 



+ 5 J]Q J rk £ J D Ja {u)I k ; k {u J ) 



(B8) 



Since these tensor elements are given in terms of an arbitrary 
basis for the donor and acceptor electronic states, we may 
write them in terms of the ICC states for which J is diago- 
nal. Considering the elements that affect populations (k = k' 
for acceptor, j — f for donor), it is then evident that in the 
ICC representation, the factors Jj ok restrict nonzero contri- 
butions from donor states j to only those deriving only from 
the coupled donor ICC state \Dj ){Dj \- Each of these terms 
also includes a sum over other inter-complex coupling matrix 
elements and lineshapes, but these only affect the magnitude 
of the allowed transitions. In the case where there is only a 
single nonzero ICC coupling, Eq. ( |B3| l thus reduces to Eq. ([5]) 
of the main text. 

Note that since our donor and acceptor lineshapes (Eq. ( |B5| ) 
and Eq. ( |A4| i) take identical forms, resulting forward and 
backward transfer rates will be equal and thus may not neces- 
sarily respect detailed balance. Therefore we do not calculate 
actual rates using Eq. 



Appendix C: Proof that classical transport is unbiased 

For simplicity, we restrict an excitation initially at the lower 
energy site (exciton) of each dimer to two possible jumps: 
with probability p to the higher energy site (exciton) of the 
same dimer, or with probability 1 — p to the higher energy 
site (exciton) of the neighboring dimer to the left. For an ex- 
citation at the higher energy site (exciton), by detailed bal- 
ance the rate of transitions is scaled by the Boltzmann factor 
e /3AB Th us although there is a longer delay before a jump, 
the probability of an inter-dimer jump is still p, and 1 — p for 
the intra-dimer jump, now to the right. Since every jump al- 
ternates between a high and low energy state, on average the 
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random walk is stationary. If non-nearest-neighbor jumps are 
allowed to states on other dimers with the same energy, then 
detailed balance requires that the probability of a jump to the 
left and the reverse jump to a right are equal, so the random 
walk is still unbiased. 



Appendix D: Ratchet methods 

a. Propagation of coherence With weak coupling be- 
tween different dimers, inter-dimer transfer should follow the 
principles of our theory of propagated coherence described 
in Section EL Here we restrict the inter-dimer coupling to be 
between nearest neighbors, to simplify the singular value de- 
composition. This is a reasonable approximation for realistic 
dipole-dipole couplings in light harvesting arrays because the 
1 /r 3 scaling ensures a rapid fall-off with inter-chromophore 
distance. Accordingly, an ICC analysis tells us that after an 
inter-dimer transfer the dynamics will be reset with the initial 
condition on the site in the dimer nearest the side from which 
the excitation was received. Thus if an excitation is received 
from the dimer to the left (right), we restart dynamics the ini- 
tial condition is on the left (right) site of the new dimer. For 
the complex consisting of the nth dimer, this corresponds to 
the explicit inter-complex coupling matrix (in the ICC basis) 



^2 J\n,e)(n-e, 

£S{-1, + 1} 



(Dl) 



where J is an arbitrary inter-dimer coupling strength and 
\n, e) is the state corresponding to occupation of the right 
(e = — 1) or left (e = +1) site of dimer n. We also assume 
that upon excitation at a site in a new dimer, the baths of the 
donating chromophore will instantaneously revert to thermal 
equilibrium. Accordingly, since the chain of dimers is peri- 
odic, we can build overall dynamics in this manner from full 
quantum calculation of just four time-dependent transfer rates 
corresponding to left or right transfer to a neighboring dimer 
from each of the two initial conditions on sites. 

b. Scaled 2CTNL calculations For computational feasi- 
bility, we based our calculations of transfer rates on scaling 
the results of 2CTNL simulations on a three dimer (six site) 
system. We may denote the left, central and right dimers 
as — 1,0, +1, respectively. We use the 2CTNL method be- 
cause it accurately models dynamics under both strong and 
weak environmental coupling |22|. Since we need to cal- 
culate rates neglecting back-transfer while these simulation 
methods describe full system dynamics, we calculate the dy- 
namics here with the inter-dimer coupling Jo set to be very 
small so that back-transfer was negligible. To begin, we 
need cumulative transition probabilities F® s (t) for the tran- 
sition from initial conditions e 6 {+1,-1} to neighboring 
dimer S 6 {+1,-1} at time t. From our simulations on 
the three dimer chain, the quantity F® s (t) is simply the to- 
tal population at time t on dimer 6 obtained by starting with 
initial condition on site e of the central dimer. The proba- 
bility density as a function of t, which is the transition rate, 
is then given by f® s (t) = §fF^ s (t) an( ^ evaluated numeri- 
cally. We then scale the transfer rate to obtain the rescaled rate 
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FIG. Al: Inter-complex transfer rates f e s(t) with initial conditions 
(a) e = +1 and (b) e = —1 for steps 6 = +1 (solid lines) and 8 = 
— 1 (dashed lines). These rates are derived from 2CTNL simulations 
for a chain of three dimers with correlation time 50 fs and no spatial 
correlations, as described in the text. 



fes(t) = {J I Jo) 2 i corresponding to coupling J. The 
rescaled cumulative transition probability is obtained by nu- 
merical integration, F c s{t) = J Q drf^ir). This scaling of 
the transfer rate assumes that to lowest order in J the transfer 
rate is proportional to J 2 , as given by Eq. ([TJ. For our param- 
eters, we found that the scaled transfer rate f e $ (t) does indeed 
converge as Jo — > and that using a value Jo = 1 cm -1 was 
sufficiently small to make any error negligible. This method 
neglects the effects of excitation loss on the dynamics of the 
donating dimer, which is reasonable to lowest order in J. We 
simulated the three dimer chain using two such 2CTNL cal- 
culations, one for each initial condition on a specific site of 
the central dimer. Calculations including spatially correlated 
baths on each dimer were performed as described previously 
|9). In principle, one could perform calculations taking into 
account static disorder, but we do not expect static disorder 
would influence our qualitative findings since the primary ef- 
fect of disorder is to limit derealization and our scaling pro- 
cedure already constrains exciton derealization to individual 
dimers. 



Figure Al gives an example of these transfer rates f e s(t). 
The transfer rates oscillate, corresponding to quantum beat- 
ings in the donor, but eventually converge to the same equi- 
librium rate, as required to satisfy detailed balance. However, 
at early times, the left (5 = +1) and right (<5 = —1) transfer 
rates are not equal, oscillating out of phase. When averaged 
over the limiting distribution ir e of the initial condition this 
gives rise to the marked short time asymmetry in the left and 
right inter-dimer transfer rates that is shown in Figure |4|b) 
of the main text. This asymmetry, although small in absolute 
terms, is amplified by being repeated over many hops between 
dimers and is responsible for the asymptotic bias of the ran- 
dom walk. 

c. Generalized random walk With each transition to a 
neighboring dimer only depending on the initial conditions at 
each dimer, the dynamics constitute a type of Markov chain 
controlled random walk known as a semi-Markov process 
ll38ll . In each step of the random walk, we start with a state 
of the form |n, e) denoting occupation of the right (e = — 1) 
or left (e = +1) site of the nth dimer. The cumulative tran- 
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sition probabilities F$ e (t) for transitioning from |n, e) to the 
left or right dimer ICC acceptor state \n + S, 5) for 5 = ±1 
(see Eq. ( |D1[ )) are determined by the rescaled 2CTNL calcu- 
lations as described above. The update ef = S is the initial 
condition for the excitation on the new dimer following from 
our ICC analysis. The walk is memory-less in terms of a two- 
dimensional clock variable (n, t) denoting "space-time" posi- 
tion, but the initial condition at each dimer nevertheless func- 
tions as an additional "coin" degree of freedom e that controls 
the likelihood of jumps to a new clock state (n', t'). 

d. Monte-Carlo algorithm Given the transition proba- 
bilities F e s (t) for this random walk, we used two techniques 
to calculate the long time behavior of the overall random walk. 
First, we performed Monte-Carlo simulations of the evolu- 
tion for a total time T by averaging over trajectories of many 
jumps. We start by setting the clock to the state (n, t) = (0, 0) 
and the coin to e = +1. We sample from the distribution of 
possible space -time shifts £ e g — (S,t e s) by choosing a pair 
(til, 1*2 ) °f independent uniformly distributed random num- 
bers between and 1. For convenience, we define the final 
transition probability p e g = lim^oo F e g(t). If ui < p e ,+i, 
we choose 6 = +1 for this jump; otherwise, S = — 1. The 
time t e g it takes for this jump is determined by numerically 
solving the equation U2 = F e g(t e g)/p e g for t. We then update 
the clock (n', t') = (n + 6,t + t e g) and the coin e' = 5. This 
process is repeated until the t + t € g > T, at which point we 
record the location of the previous dimer n as the final state of 
that trajectory. The probability density of the final distribution 
over dimers is derived by binning over many such trajectories 
(^5000). Empirically, our Monte-Carlo simulations suggest 
that the distribution converges to a normal distribution char- 
acterized by its mean and variance, as expected from a central 
limit theorem for weakly dependent variables [39|. 

e. Analytic model Second, we calculated the mean and 
variance of final distribution analytically in the asymptotic 
limit of the total walk time T — > oo, using the method 
suggested in Ref. [34] Since the results of these calcula- 
tions agreed with the Monte-Carlo simulations but were much 
faster, we use this second method for the plots in Figure|4] To 
begin, we calculate the moments of the transition time t e g for 
each jump e — > 5, 



1 [°° 
E(t 6 «) = — / tf eS (t)dt 

PeS Jo 

i r°° 

E(&) = — / t 2 fMt)dt 
PtS Jo 



(D2) 
(D3) 



by numerical integration. Now, note that transitions between 
coin states can be described as a Markov chain with transition 
matrix P with entries given by the final transition probabilities 
p e s = lim^oo F e s(t) as defined above, i.e., 



P = 



P+i,+i P+i-i 
P-i,+i P-i-i 



(D4) 



Accordingly, the limiting distribution 7r e over the coin space 
is given by the left eigenvector of P with eigenvalue 1, that is, 
the solution ir of the equation 7r e < = Yle ^^P^' ■ The quantity 
Ti € p e s gives the probability of the step e — > 5 in the limiting 



distribution. Recalling the definition of the space-time shift 
its = (S, t e g) associated with the step e — > 6, we obtain an 
average space-time shift £ £ g — (5, ~E(t e s)) for this step. Since 
the coin will converge to the limiting distribution 7r e , we now 
obtain the average space-time shift over all steps as 



(D5) 



e<5 



Now let riT denote the spatial position of the random walk 
after a total time T. This random walk is the sum of T/i 
independent steps on average, each of which has an average 
spatial shift n. Since the expectation adds linearly, we then 
obtain the average position of the overall walk as 



E(n T ) = nT/t. 



(D6) 



Figure |4]plots the corresponding drift velocity v = E(n^)/r. 

Given that our random walk appears to converge to a nor- 
mal distribution, we can fully characterize the distribution 
with its mean, calculated above, and its variance. As a practi- 
cal matter, the variance indicates the width of the distribution 
and thus determines whether or not a non-zero drift velocity 
would be observable experimentally. To calculate the vari- 
ance, we consider two sources of space-time deviations, 



VtS = £e<5 — £e<5 



(D7) 
(D8) 



corresponding to deviations ri e s of the space-time shift of a 
particular transition from its average value, and deviations fi e g 
of the average space-time shift for a particular transition from 
the average space-time shift over all transitions. Since suc- 
cessive steps are weakly correlated by the ICC conditions, the 
latter quantity must be averaged over all possible steps in all 
possible trajectories. We therefore define p, as the single step 
average obtained by summing fj, e g over all possible sequential 
steps: 



(i = lim -y^UaSi 

i=l 



(D9) 



Note that for a standard Markov chain with no correlation be- 
tween steps, this single step average reduces to \i. Assuming 
the covariation between r\ t g and p, is small, we can combine 
them to calculate the overall space-time covariance matrix, 



Var(£) = Var(£ - £) 
= Var(?7 + p) 
= Var(r;) + Var(/l). 



(D10) 
(Dll) 
(D12) 



Averaging over the limiting distribution 7r c , we find for the 
first contribution to the variance, 



Var(?y) = 7r £ p e5 Var(r? e5 ) 



eS 



where 



Var(77 e<5 ) 





Var(^) 



(D13) 



(D14) 
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with Var(t e 5) = E(i^ (5 ) — E(i e a) 2 . Now consider the second 
contribution to the variance, Var(/2). The variance of the sin- 
gle step average ft,, Eq. ( |D9| l, introduces a double sum over 
products of deviations fi e g. With a correlated random walk, 
the products of deviations at different space-time values are 
also correlated and hence evaluation of these requires enumer- 
ation of all possible jumps connecting them, where these are 
determined by the transition probability matrix P, Eq. ( |D4| i. 
This enumeration, which constitutes a multi-state generaliza- 
tion of the variance for weakly dependent processes ll39l . is 
given explicitly by 



Var(^) = ^7T, 



+ ^ePeSPs^ Ppa^sPp* + f^Ves], (D15) 

eSptr m>0 

where pf^ — (P m )s P and we sum each variable e,S,p,a 
over ±1. The second term in Eq. D15 sums up all contri- 
butions that m-steps apart, where these are specified by the 
Chapman-Kolmogorov equation ETI . We note that for con- 
venience, instead of explicitly performing the sum over m, 
one can equivalently replace the term J2 m vfp* m e 1 ua " 
tion above with Qs p B4l . where Q = (1 — P*)^ 1 and 
P* is the non-equilibrium portion of P, that is, with entries 
p* s — p eS — 7T C . Combining Eqs. ( |D 1 3| > and ( |D15| l into ( |D 1 3| > 
yields a space-time covariance matrix Var(£) for the two di- 
mensional shift variable £. This covariance matrix has explicit 
entries, 



Vax(0 



/ Var(n) Cov(n, t) 
\Cov(n, t) Var(i) 



(D16) 



To calculate the final spatial variance Var(nr), we must take 
into account the uncertainty associated with the number m of 
discrete hops that happened in time T, in addition to the un- 
certainty over n. To correctly incorporate both contributions, 
we calculate the variance of the spatial displacement n over a 
single coin shift over the full two-dimensional coin space, 

Var( n> t) (n) = Vax(„ )t ) (n - nt/t) 

= Var(n) - 2(n/t) Cov(n, t) + (n/t) 2 Var(i) 
= (l > -n/t)Var(0(l,-n/t) T , (D17) 

where in the first step we subtracted the average value of n 
over the coin space. We write the variances over the full coin 
space (n, t) to emphasize that they are distinct from terms like 
Var(n), which is only over the spatial degree of freedom n. 
Since the variance adds linearly over m « T/t independent 
steps, we obtain the variance of the distribution after time T 
as 

Var(n T ) = (1,-n/t) Var(£)(l, -n/£) T j. (D18) 

The diffusion coefficient for the walk is then given as D = 
Var(n r )/2T. 

Figure |A2| plots the full results of scans over correlation 
time and cross-correlation coefficients used to create Figure [4] 



of the main text. We see that the width of the excitation trans- 
fer distribution is approximately constant over all parameter 
choices at about 60 rim after 1 ns, and that the asymmetry be- 
tween initial conditions Att = tt + i — 7r_i accounts for most 
of the variation in drift velocity. Figure Qb) of the main text 
is a plot of the relative transfer rate asymmetry for the limiting 
distribution of the initial condition 7r e , 



A(t) = 



(D19) 



where the sums are over e, S £ { — 1, +1} as usual. 



Appendix E: FMO Hamiltonian and singular value 
decompositions 



The FMO complex exists in a trimer arrangement, where 
each monomer contains 7 bacteriochlorophyll molecules, and 
three additional BChl molecules (termed the 8th Bchl for each 
of the three monomers) are each located between a distinct 
pair of monomers ll24l . In this paper, we use a Hamilto- 
nian for a monomer of the FMO complex of C. tepidum cal- 
culated by Adolphs and Renger lfT4l . augmented by dipole- 
dipole couplings to the 8th BChl site calculated using their 
same methodology with the crystal structure of Tronrud et al. 
Il24ll . We assign each of the three BChl 8 pigments to the 
monomers with which they have the strongest dipole-dipole 
coupling. The Hamiltonian matrix is given below in units of 
cm -1 above 12 210 cm -1 , where elements of the matrix are 
indexed according to site from 1 to 8: 
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The energy of site 8 is marked with a question mark to indicate 
that it is unknown, since it has not been calculated. Accord- 
ingly, our simulations of the full FMO complex use only the 
portion of this Hamiltonian for sites 1-7, as in previous studies 

El 

To determine donor and acceptor ICC states for a given cou- 
pling matrix J, we perform the singular value decomposition 
J = U D 3ll\ = ^2iJi\Di)(Ai\ as described in Section [n] 
Here are the results of two examples we use with our FMO 
Hamiltonian. Let the notation denote the coupling ma- 
trix from the donor (D) rows and the acceptor (A) columns of 
Eq. (JETJ. As plotted in Figure [T] for the coupling from site 8 
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FIG. A2: Full results of simulations used for analysis of the unidi- 
rectional random walk. The parameters of each dimer match that 
of sites 1-2 in FMO, as described in the main text. We used a De- 
bye spectral density with a reorganization energy of 35 cm" 1 and 
variable bath and spatial correlations, as indicated on the figure. 
For the left panels we have variable time correlation and no spa- 
tial correlations. For the right panels we vary the spatial correla- 
tion (see Ref. 9 1 and fix the correlation time at 50 fs. ( a,b) Drift 
velocity, from Eq. ( |D6| >. (c,d) Standard deviation a = ^/Var(riT) 
of the walk at 1 ns, from Eq. ( |D18} . (e,f) Coherence time r, from 
a least squares fit of the exponential decay of excitonic coherence, 
|p eie J ~ Ae~ t/T + B: we evaluate this here and in Figure 0for 
t > 100 fs to exclude non-exponential decay. (g,h) Asymptotic 
transfer asymmetry A7r = n+i — 7r_i indicating the overall pref- 
erence for right over left transfer. 



to sites 1-7, we have J^" 7 = J*\D*)(A*\ with 



0.912" 
0.158 
0.031 

\A*) = 0.043 , (E2) 
-0.105 
0.229 
0.275 



J* = 41.9 ID* 



with entries (i\A*) for states \i) = |1), . . . , |7). With only 
a single donor site, the ICC acceptor (donor) state from the 
singular value decomposition is as simple as the normalized 
vector corresponding to the dipole-dipole matrix. Since occu- 
pation probabilities correspond to these amplitudes squared, 
the acceptor among site 1-7 is mostly (83%) on site 1. 

A less trivial example is given by considering the 1-2 dimer 
as a donor to and acceptor from the remainder of the complex 
3-8. In this case, we have Jf if = J2i=i 2 Ji\Di){Ai\ with 



Ji = 43.6 \D X ) 



J 2 = 34.3 \Do) = 



-0.2971 
0.085 

-0.969] 1 4 \ -0.153 

-0.247 1 x ' ~ 0.238 
0.196 
-0.887. 

-0.832 
-0.274 

0.247 1 , r , 0.028 
-0.969 1 2 > ~ -0.432 
-0.193 
0.089 



(E3) 



(E4) 



with entries corresponding to states \i) in ascending order. 
Thus the coupling in the ICC basis is mostly from site 8 (78%) 
to site 1 (94%), and from site 2 (94%) to site 3 (69%). If we 
omit site 8 from the acceptor, these acceptor and donor states 
are modified as follows: 



Ji = 19.7 |Di 



J 2 = 34.4 |D 2 ) = 



0.995 
0.099 



0.099 
-0.995 



\Ao) = 



0.433 • 
-0.257 

0.342 
-0.633 
-0.479. 

-0.876" 
-0.254 
-0.001 
-0.381 
-0.153 



(E5) 



(E6) 



We use these ICC states in Fig. [5] since only sites 1-7 are in- 
cluded in the 2CTNL simulation, as the site energy of the 8th 
BChl is unknown (see Eq. (JET}), and as it is furthermore un- 
clear whether this 8th BChl is present in all cases in the natural 
system 11241 . 
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